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Abstract 

We present a systematic theory of Coulomb-induced correlation effects in the nonlinear optical 
processes within the strong-coupling regime. In this paper we shall set a dynamics controlled 
truncation scheme [lj] microscopic treatment of nonlinear parametric processes in SMCs including 
the electromagnetic field quantization. It represents the starting point for the microscopic approach 
to quantum optics experiments in the strong coupling regime without any assumption on the 
quantum statistics of electronic excitations (excitons) involved. We exploit a previous technique, 
used in the semiclassical context, which, once applied to four-wave mixing in quantum wells, allowed 
to understand a wide range of observed phenomena [2fl. We end up with dynamical equations 
for exciton and photon operators which extend the usual semiclassical description of Coulomb 
interaction effects, in terms of a mean-field term plus a genuine non-instantaneous four-particle 
correlation, to quantum optical effects. 
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I. INTRODUCTION 

Since the early Seventies |3| researchers have been exploring the possible realization of 
semiconductor-based heterostructures, devised according to the principles of quantum me- 
chanics. The development of sophisticated growth techniques started a revolution in semi- 
conductor physics, determined by the possibility of confining electrons in practical structures. 
In addition, the increasing ability in controlling fabrication processes has enabled the ma- 
nipulation of the interaction between light and semiconductors by engineering, in addition 
to the electronic wave functions, the light modes. 

Entanglement is one of the key features of quantum information and communication tech- 
nology [4| and a hot topic in quantum optics too. Parametric down-conversion is the most 
frequently used method to generate highly entangled pairs of photons for quantum-optics ap- 
plications, such as quantum cryptography and quantum teleportation. Rapid development 
in the field of quantum information requires monolithic, compact sources of nonclassical 
photon states enabling efficient coupling into optical fibres and possibly electrical injec- 
tion. Semiconductor-based sources of entangled photons would therefore be advantageous 
for practical quantum technologies. The strong light-matter interaction in these systems 
gives rise to cavity polaritons which are hybrid quasiparticles consisting of a superposition 
of cavity photons and quantum well (QW) excitons 5j. Demonstrations of parametric am- 
plification and parametric emission in semiconductor microcavities (SMCs) with embedded 
QWs , Q], together with the possibility of ultrafast optical manipulation and ease of 
integration of these microdevices, have increased the interest on the possible realization of 
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nonclassical cavity-polariton states 
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[12], |l3j]. In 2004 squeezed light generation in 
14[ | . In 2005 an experiment 



SMCs in the strong coupling regime has been demonstrated 
probing quantum correlations of (parametrically emitted) cavity polaritons by exploiting 
quantum complementarity has been proposed and realized jllj ]. Specifically, it has been 
shown that polaritons in two distinct idler modes interfere if and only if they share the 
same signal mode so that which-way information cannot be gathered, according to Bohr's 
quantum complementarity principle. 

Laser spectroscopy in semiconductors and in semiconductor quantum structures has been 
greatly used because exciting with ultrashort optical pulses in general results in the creation 
of coherent superpositions of many-particle states. Thus it constitutes a very promising 
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powerful tool for the study of correlation and an ideal arena for semiconductor cavity quan- 
tum electrodynamics (cavity QED) experiments as well as coherent control, manipulation, 
creation and measurement of non-classical states 0,0 , ]l5|, ]l6| . The analysis of nonclassical 
correlations in semiconductors constitutes a challenging problem, where the physics of inter- 
acting electrons must be added to quantum optics and should include properly the effects of 
noise and dephasing induced by the electron-phonon interaction and the other environment 



channels 



171 ] . The nonlinear optical properties of exciton-cavity system play a key role in 



driving the quantum correlations and the nonclassical optical phenomena. The crucial role 
of many-particle Coulomb correlations in semiconductors marks a profound difference from 
the nonlinear optics of dilute atomic systems, where the optical response is well described 
by independent transitions between atomic levels, and the nonlinear dynamics is governed 
only by saturation effects mainly due to the balance of populations between different levels. 

The Dynamics Controlled Truncation Scheme (DCTS) provides a (widely adopted) start- 
ing point for the microscopic theory of the light-matter interaction effects beyond mean-field 
151 ] . supplying a consistent and precise way to stop the infinite hierarchy of higher-order 
correlations which always appears in the microscopic approaches of many-body interacting 
systems without need to resort to any assumption on the quantum statistics of the quasi- 
particle arising in due course. By exploting this scheme, it was possible to express nonlineari- 
ties originating from the Coulomb interaction as an instantaneous mean-field exciton-exciton 
interaction plus a noninstantaneous term where four-particle correlation effects beyond men- 
field are contained entirely in a retarded memory function In 1996 the DCTS was ex- 
tended in order to include in the description the quantization of the electromagnetic field 
and polariton effects [3]. This extension has been applied to the study of quantum opti- 
cal phenomena in semiconductors and it was exploited to predict polariton entanglement 
[13I . The obtained equations showed that quantum optical correlations (as nonlinear optical 
effects) arise from both saturation effects (phase-space filling) and Coulomb induced correla- 
tions due to four-particle states (including both bound and unbound biexciton states). The 
dynamical equations included explicitly biexciton states. The structure of those equations 
didn't allow the useful separation of Coulomb interaction in terms of a mean-field interaction 
term plus a noninstantaneous correlation term performed in the semiclassical description. 

In this paper we shall set a DCTS microscopic treatment of nonlinear parametric pro- 
cesses in SMCs including the light-field quantization. It represents the starting point for the 



3 



microscopic approach to quantum optics experiments in the strong coupling regime. For this 
purpose we shall exploit a previous technique [2j which, once applied to four-wave mixing 
in QWs, it allowed to understand a wide range of observed phenomena. Here all the ingre- 
dients contributing to the dynamics are introduced and commented. We shall give in great 
details the manipulations required in order to provide an effective description of the non- 
linear parametric contributions beyond mean-field in an exciton-exciton correlation fashion. 
In particular we derive the coupled equations of motion for the excitonic polarization and 
the intracavity field. It shows a close analogy to the corresponding equation describing the 
semiclassical (quantized electron system, classical light field) coherent response in a QW 
j^], the main difference being that here the (intracavity) light field is regarded not as a driv- 



ing; external source but dynamical field 



19| . This correspondence is a consequence of 



the linearization of quantum fluctuations in the nonlinear source term here adopted, na mely 
the standard linearization procedure of quantum correlations adopted for large systems [20| . 
However the present approach includes the light field quantization and can thus be applied 
to the description of quantum optical phenomena. Indeed, striking differences between the 
semiclassical and the full quantum descriptions emerge when considering expectation val- 
ues of exciton and photon numbers or even higher order correlators, key quantities for the 
investigation of coherence properties of quantum light 11]. This is the main motivation 
for the derivation of fully operatorial dynamical equations, within such lowest order non- 
linear coherent response, we address in the last section. The results here presented provide 
a microscopic theoretical starting point for the description of quantum optical effects in 
interacting electron systems with the great accuracy accomplished for the descri ptio n of 
the nonlinear optical response in such many-body systems, see e.g. 



references therein. The proper inclusion of the detrimental environmental interaction, an 
important and compelling issue, is left for a detailed analysis in another paper of ours 23]. 

In Section [III the generality of the coupled system taken into account are exposed, here 
all the ingredients contributing to the dynamics are introduced and commented. The linear 
and the lowest nonlinear dynamics is the subject of Sec. IIII1 whereas in Sec. IIVI we shall 
give in great details the manipulations required in order to provide an effective description of 
the nonlinear parametric contributions beyond mean-field in an exciton-exciton correlation 
fashion. In Sec. |V] the operatorial equations of motion for exciton and intracavity photon 
operators are derived. 



4 



II. THE COUPLED SYSTEM 



The system we have in mind is a semiconductor QW grown inside a semiconductor planar 
Fabry-Perot resonator. In the following we consider a zinc-blende-like semiconductor band 
structure. The valence band is made from p-like (I = 1) orbital states which, after spin- 
orbit coupling, give rise to j = 3/2 and j = 1/2 decoupled states. In materials like GaAs, 
the upper valence band is fourfold degenerate (j = 3/2), whereas in GaAs-based QWs the 
valence subbands with j = 3/2 are energy splitted into two- fold degenerate heavy valence 
subbands with j z = ±3/2 and light lower energy subbands with j z = ±1/2. The conduction 
band, arising from an s-like orbital state (1=0), gives rise to j = 1/2 twofold states. In 
the following we will consider for the sake of simplicity only twofold states from the upper 
valence and lowest conduction subbands. As a consequence electrons in a conduction band 
as well as holes have an additional spin-like degree of freedom as electrons in free space. 
When necessary both heavy and light hole valence bands or subbands can be included in the 
present semiconductor model. Only electron-hole (eh) pairs with total projection of angular 
momentum a = ±1 are dipole active in optical interband transitions. In GaAs QWs photons 
with circular polarizations a = — (+) excite electrons with j e z = ±1/2 (j'J = —1/2) and holes 
with j h z = -3/2 (j h z = 3/2). We label optically active eh pairs with the same polarization 
label of light generating them; e.g. o = ±1 indicates an eh pair with j'J = —1/2 and 
3 h z = 3/2. 



We start from the usual model for the electronic Hamiltonian of semiconductors |15l. \24\. 
It is obtained from the many-body Hamiltonian of the interacting electron system in a lattice, 
keeping explicitly only those terms in the Coulomb interaction preserving the number of 
electrons in a given band, see Appendix [Bl The system Hamiltonian can be rewritten as 

H e = H + Vboul = J2 En *\ Na ) ( Na I > W 

Na 

where the eigenstates of H e , with energies E^a = fajjNa, have been labelled according to the 
number N of eh pairs. The state | iV = 0) is the electronic ground state, the N = 1 subspace 
is the exciton subspace with the additional collective quantum number a denoting the exciton 
energy level n, the in-plane wave vector k and the spin index a. When needed we will adopt 
the following notation: a = (n, k) with k = (k, a). In QWs, light and heavy holes in valence 
band are split off in energy. Assuming that this splitting is much larger than kinetic energies 
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of all the involved particles and, as well, much larger than the interaction between them, 
we shall consider only heavy hole states as occupied. On the contrary to the bulk case, in a 
QW single particle states experience confinement along the growth direction and subbands 
appear, anyway in the other two orthogonal directions translational invariance is preserved 
and the in-plane exciton wave vector remains a good quantum number. Typically, the energy 
difference between the lowest QW subband level and the first excited one is larger than the 
Coulomb interaction between particles, and we will consider excitonic states arising from 
electrons and heavy holes in the lowest subbands. 

Eigenstates of the model Hamiltonian with N=l (called excitons) can be created from 
the ground state by applying the exciton creation operator: 

\lnak) = Bl k \N = 0) , (2) 

which can be written in terms of electrons and holes operators as 

ALrk = < ^™crk'cl-,k'+7 ?e k/2^!-,-k'+r; h k/2 ' (3) 

k' 

here $^ CTk / is the exciton wave function, being k the total wave vector k = k e + k^, and 
k' = r] e k e — rjhkh with 77( e ,7i) = wi( e ,fe) /( m (e) + m (h)) {jn e and nih are the electron and hole 
effective masses). These exciton eigenstates can be obtained by requiring the general one eh 
pair states to be eigenstates of H e : 

iTJlncrk) = huin^lnak} , (4) 



and projecting this secular equation onto the set of product (eh) states \k e , kh) = c k e ^k h I 0)> 
(see Appendix IB1 for details): 



y~]{(K> K I He I h, h) ~ ftunakh> e k> h ,k e k h )(k e , k h \lnak) = . (5) 

Thus, having expressed the correlated exciton state as a superposition of uncorrelated prod- 
uct states, 



\lnak) = ^ I (k e ,k h \lnak) \\k e ,k h ) 

k e ,k h \ / 



(6) 



the scalar products, coefficients of this expansion, represent nothing but the envelope func- 
tion $j£ ak , of the excitonic aggregate being the solution of the corresponding Schrodinger 
equation (jSJ). It describes the correlated eh relative motion in k-space. In order to simplify 
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a bit the notation, the spin convention in Eq. ([3]) has been changed by using the same label 
for the exciton spin quantum number and for the spin projections of the electron and hole 
states forming the exciton. 

The next relevant subspace (N = 2) is the biexciton one, spanning all the states with 2 
eh pairs. It seems worth noting that the above description of eh complexes arises from the 
properties of quantum states and, once fixed the system Hamiltonian, no approximations 
have been introduced insofar. Indeed such a property hold for any N eh pairs aggregate and 
we will give a full account of it in Appendix [B] 

The eigenstates of the Hamiltonian H c of the cavity modes can be written as | n, A) where 
n stands for the total number of photons in the state and A = (k^ai; k„, a n ) specifies 
wave vector and polarization a of each photon. Here we shall neglect the longitudinal- 
transverse splitting of polaritons 25| originating mainly from the corresponding splitting of 
cavity modes. It is more relevant at quite high in-plane wave vectors and often it results to 
be smaller than the polariton linewidths. The present description can be easily extended to 
include it. We shall treat the cavity field in the quasi- mode approximation, that is to say 
we shall quantize the field as the mirrors were perfect and subsequently we shall couple the 
cavity with a statistical reservoir of a continuum of external modes. This coupling is able 
to provide the cavity losses as well as the feeding of the coherent external impinging pump 
beam. The cavity mode Hamiltonian, thus, reads 

H c = Yl > ( 7 ) 

k 

with the operator a\ which creates a photon state with energy hw c k = h(cu 2 xc + v 2 \k.\ 2 Y^ 2 , 
v being the velocity of light inside the cavity and k = (a, k). The coupling between the 
electron system and the cavity modes is given in the usual rotating wave approximation 



is, 



26] 



H i = ~ E V nA B nk + H.c. , (8) 
nk I 1 

V n k is the photon-exciton coupling coefficient enhanced by the presence of the cavity [27J 
set as V nt k = ^a-V^40n(j( x = 0), the latter being the real-space exciton envelope function 
calculated in the origin whereas A is the in-plane quantization surface, V a is proportional to 
the interband dipole matrix element. Modeling the loss through the cavity mirrors within 
the quasi-mode picture means we are dealing with an ensemble of external modes, generally 
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without a particular phase relation among themselves. An input light beam impinging on 
one of the two cavity mirrors is an external field as well and it must belong to the family of 
modes of the corresponding side (i.e. left or right). Being coherent, it will be the non zero 



expectation value of the ensemble. It can be shown 



18 



231 ] that for a coherent in put beam, 



the driving of the cavity modes may be described by the model Hamiltonian 18|, [23] 



H p = it c ^2(E k al - ££a k ) 



(9) 



where (E£) is a C-number describing the positive (negative) frequency part of the co- 
herent input light field amplitude. 



III. LINEAR AND NONLINEAR DYNAMICS 

The idea is not to use a density matrix approach, but to derive directly expectation values 
of all the quantities at play. The dynamics is described by "transition" operators (known as 
generalized Hubbard operators): 

X N , a . M ,f3=\ N,a){M,/3 | 
Yn,\;m,» =1 n, \) (m, fx I . (10) 

The fundamental point in the whole analysis is that, thanks to the form of the interaction 
Hamiltonian Hj and thanks to the quasiparticle conservation the free Hamiltonians possess, 
we can use the so-called dynamics controlled truncation scheme, stating that we are facing a 
rather special model where the correlation have their origin only in the action of the electro- 



magnetic field and thus the general theorem due to Axt and Stahl [1] holds. For our purpose 
we will need its generalization in order to include the quantization of the electromagnetic 
field [l8[], it reads: 

io 



(N+M+n+m+2i) 



i=0 



+0(-E ,(Ar+M+n+m+2io+2 )) (11) 



i.e. the expectation value of a zero to N-pair transition is at least of order N in the 
external electromagnetic field. There are only even powers because of the spatial inversion 
symmetry which is present. Once a perturbative order in the external coherent fields is 
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chosen, Eq. (fTTj) limits the expectation values to take into account, thus providing a precise 
way to truncate the hierarchy of equations of motions. 
The exciton and photon operators can be expressed as 

k;(n k +l)k 

n>l 

B nk = X ;lnk+ (Na \ B n k I (N + l)(3) X Na , {N+1)f3 , (12) 

N>l,a/3 

where in writing the photon expansion we omitted all the states not belonging to the fe-th 



mode which add up giving the identity in every Fock sector 28]. 

The equation of motion for the generic quantity of interest Xjsr^M^Yn^im,^ is reported in 
Appendix [A] In the Heisenberg picture we start considering the equation of motion for the 
photon and exciton operators, once taken the expectation values we exploit theorem (TTTj) 
retaining only the linear terms. With the help of the generalized Hubbard opertors all this 
procedure may be done by inspection. The linear dynamics for ( a k = ( Y 0; x nk )^ and 
( B nk = ( X . lnk reads: 

| ( a, )« = -ifi>° ( a k >« + iJ2 V ~f < B nk + t c f (13) 

n 

j t ( B nk = -iO lnk ( B nk + i^(a k . (14) 

In these equations u) k = uj k — ij k , where ^ k is the cavity damping, analogously u)\ nk = 
^lnk ~ *r x and u>2/3 — ^2/3 — ff xx . The dynamics up to the third order is a little bit more 
complex, we shall make extensively use of flAlj) (in the following the suffix + ( n ) stands for 
"up to" n-th order terms in the external electromagnetic exciting field). With Eq. (fi~2l) the 
exciton and the photon expectation values can be expanded as follows: 

( B nk ) +(3) = ( X , lnk ) +(3) + ^(la | B nk | 2(3) ( X la , w ) (3) , (15) 

a/3 

( a k >+ (3) = ( %, lk > +(3) + V2 ( Y lk , 2k ) (3) . (16) 
With a bit of algebra we obtain 

| ( a k > + < 3 > = -iut ( a k > + < 3 > +iJ2 Y f< B nk >+< 3 > + , (17) 

n 
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^ / D \+(3) ■- / 6 \+(3) i -Vnk I - \+(3) . 

-J t \ B nk) = -^lnfc ( B nk ) y ' + I— (a k ) K ' + 



rik 



^ ^ V n r h ,(lhk I [5^,5^,] -5(^)5^) | loc)(X liik . la Y . lk , ) (3) 

n'k',a 



-i 5^(w 2j 9 - ^ifife - wir»fc)(lnfe I 5 n fc I 2/3)(l 1 ^ ;2/3 i > o;o ) (3) , (18) 

in analogy with the eqs [l^] (see also Ref. [2I). The resulting equation of motion for the 
lowest order biexciton amplitude is 

+1 Yl V nlk ,(2p\Bl k ,\ln"k")(X^ ln ,,«,% tlk ,f ) . (19) 

n'k';n"k" 

IV. COHERENT RESPONSE 

Thanks to the fact we want to treat coherent optical processes it is possible to ma- 
nipulate further the parametric contributions under two assumptions. We are addressing 
a coherent optical response, thus we may consider that a coherent pumping mainly gen- 
erates coherent nonlinear processes, as a consequence the dominant contribution of the 
biexciton sector on the third-order nonlinear response can be calculated considering the 
system quantum state as a pure state, which means the nonlinear term is regarded as orig- 
inating mainly from coherent contributions. Moreover nonclassical correlations are taken 
into account up to the lowest order. The first assumption results in the factorizations 
(X lh - k ^Y ;of } * (^ ; o> (1) (^o ;2/3 > (2) and (X m ~ k . ip Y . >lk , ) (3) * (X mk , )^ (X ;ip%,ik> f ] ■ 
The second implies (^o;i,s^o;ifc' )^ — (Xo-,ip)^ {Y .i k > } , m the nonlinear source term, 
namely the standard linearization procedure of quantum correlations adopted for large sys- 
tems [20|. Of course these two approximations can be avoided at the cost of enlarging the 
set of coupled equations in order to include the equation of motions for the resulting cor- 
relation functions. It neglects higher order quantum optical correlation effects between the 
electron system and the cavity modes leading to a renormalization of the biexciton dynamics 
with intriguing physical perspectives. However for extended systems, like QWs in planar 
microcavities, these are effects in most cases of negligible impact, on the contrary in fully 
confined geometries such as cavity embedded quantum dots they could give significant con- 
tributions. In the end, within such a coherent limit, we are able to describe the biexciton 
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contribution effectively as an exciton-exciton correlation [2]. The resulting equations for 
the coupled exciton an cavity-field expectation values coincide with those obtained within 
a semiclassical theory (quantized electron-system and classical cavity field). Nevertheless 



completely different results can be obtainec 



or for higher order correlation function 13, 



or exciton or photon number expectation values 



261 ] . In the next section we will derive operator 



equations useful for the calculation of such correlation functions. After the two approxima- 
tions described above (linearization of quantum fluctuations and coherent limit), Eqs ffT8j) 
becomes 

± ( B nk >+< 3 > = -iu lnk ( B nk >+< 3 > + ^ < a k - 1 £<4*>* ( X 2 fe U ' ( 2 °) 



where 



with 



dt "" " H h 

nk 



nk;nk nk;nk nk;nk v ' 
n'k',n"k" 

n COUL(2) _ V~> (1) ,j> v( 2 ) / 91 x 

P 



C Snk" k " = V n'k>(lnk | 5 {n , k , );{nk) - [B nk ,Bl k ,\ I ln"k") (22) 



c 2!.«S.fl = ^( W2 /9 ~ ~ Wi^fcXlnfe | S nfe | 213) . (23) 



(1) 

"nk;hk;/3 

This equation is analogous to the corresponding equation describing the semiclassical (quan- 
tized electron system, classical light field) coherent x response in a QW [2j, the main 
difference being that here the (intracavity) light field is regarded not as a driving external 
source but as a dynamical field ljj. This close correspondence for the dynamics of ex- 
pectation values of the exciton operators, is a consequence of the linearization of quantum 
fluctuations. However the present approach includes the light field quantization and can 
thus be applied to the description of quantum optical phenomena. By explicit calculation 
it is easy to see that the first term in Eq. (1221) is zero unless all the involved polarization 
labels a coincide. In order to manipulate the last term we follow the procedure of Ref. ^ 
which succeeded in reformulating the nonlinear term coming from the Coulomb interaction 
as an exciton-exciton (X-X) mean-field contribution plus a correlation term driven by a two- 
exciton correlation function. Even if we are about to perform more or less the same steps 

11 



of Ref. Q| we shall provide a detailed account of all the key points of the present deriva- 
tion. A clear comprehension of these details will be essential for the extension to operatorial 
dynamical equations of the next section. 

In performing this we shall need the two identities: 

C nl,nk;P = ~ U ^' k ~ W lnfc)(!^ I B nk | 2(3) = 



H 

K(lnk | B nk (-^ - oo mk - uo lnk ) | 2(3) (24) 



and 



d 
di 



(4*) (1) (ii«) (1) e^j = (25) 

= +^(v n , v (av)W(B n ,w)W + W<ak») (1) <4w))e~ in(t ~ t ' ) , 
or 

1 d sr~^ I ,~ , m , " ,m son *'\\ i 



E (4*> (1) (^»> (1) e- ifi(t - t ' ) )=+^ E WM (1) <W> (1) e- m( ^\ 

n'k';n"k" \ ' n'k';n"k" 

(26) 

where f2 = wi n /fc/ + Wi n "fc" — 2iT x . Employing the forma solution of the biexciton amplitude 
Eq. (I19p we have: 



E c iW*°^ (2) = n T,( m ~ k i - - ^ fc ) 1 2/3) ■ (27) 

. Y^ {2( 3\Bl k ,\ln"k") f ^^^^(ovJ^CO^O^CO- 

n'k';n"k" 

We observe that the matrix elements entering the nonlinear source terms are largely inde- 
pendent on the wave vectors for the range of wave vectors of interest in the optical response. 
Neglecting such dependence we can thus exploit the identity ( 1261) . obtaining 



E (2/5 I ^ I ln " k ")\ 4j (t')(B n ^ (Oe-^(*-*')V^(«') = 

n'k';n"k" \ / 

h E /' I - "if* - uink)e-^-ti)Bl iki | i„^> e - r -C*-0 

~, nun */ — oo 



n'k';n"k" J ~°° 

1 d 

2 eft 
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where in the last lines we have resummed all the biexciton subspace by virtue of its com- 
pleteness. By performing an integration by part, Eq. fl28|) can be rewritten as 



1 



n'k';n"k" 



3 ! ( w l»'i'+ w ln"t"- 2ir x+ ir xx)( i - ! ') 



(Ink | B nk (^ - u m% - u; lnk )e-^~^Bl k ,\ ln"k") 



(B n >k>) (1) (t')(B n „ k „) (1) (t')e-^-^ 



dt\B n , k ,) {1) (t')(B n ,, k n f\t') e - iQ{ - 1 -^ — ) ■ 



df 



H, 



(Ink | B nk {-^ - u mk - uj Xnk )e-^-^B\ lk , \ ln"k") 



(29) 



^hJ2 \(lnk\B nk ( 

n'k';n"k" I 



He 

h 



"ink ~ Umk)Bl k ,\ ln"k")(B n , k , (t)(B n//k „ ) w (t) 



,(1) 



(Ink | B nk (-^ - u lfik - u^e-^-^B-lu \ ln"k") \ . 



(30) 



The first and the second term can be expressed in terms of a double commutator structure: 



B 

(Ink | B nk (-^ - u mk - uj lnk ) = (0 | [B nh [B nk , H c ]] = (0 | D hk nk , 



(31) 



where a force operator D is defined |2j and 



jL) (3 *K„'fc'+^ln''fc''- 2ir x + ir x x)(t-t') 

dt>^ 



(Ink | B nk (^ - - uj lnk )e-^-^Bl k/ \ ln"k") 



= ^|<0 I D^^^ ilk ,Bl kll | 0>e^- -«r^r-)(^) 
H / H \ 

= (° I D^^e-^^^ii^ - w ln / fc , - w ln /, fc /, - i(r xx - 2r x )J 

n n'k ,JD n"k" I U / e — 

= e »K„'fc'+ w ln''fc''- 2 ff X +irxx)(t-t')^ n " fe ''V fc '(i _ f>) _|_ 

hk,nk * ' 

+(r xx - 2r x )e i ^'^' + ^"^"- 2irx+ir -)( t - t ' ) (o | D^ nfc (t - O^v^fc" I °) > 



(32) 
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where the memory kernel reads 

<£ BV (* - = <° I - O^w* I 0) • (33) 

The usual time dependence in the Heisenberg picture is given by D(r) = e l< y H d h ) T £) e - l ( H c/ h ) T . 
Altogether, the nonlinear term originating from Coulomb interaction can be written as 

n COUL(2) _ (1) (y v(2) _ 



n'k';n"k" I 



/ ^J^C* - Ofc (0 (t')e-^W - (34) 

>/ — oo I 

(r„ - 2r x )]T f dt'<o | 4^(t - t)Bl„& n „r | o)<J3„v (t'X^W (f ) • 



n'k' 
n"k" 



For later purpose we are interested in the optical response dominate by the IS exciton sector, 
with T xx ~ 2r x in the cases of counter- and co-circularly polarized waves. Specifying to this 
case the Coulomb-induced term with Eq. (|34"|) becomes 



^<B ± +(3, -^£<4i;>* (I, C."!f , = < 35 ) 



h 

COUL 



-iO k ( B ±k -l£ W; k ' +k "^ x (5 ± k)* (1) (i)(^) (1) W(5 ±k ,) (1) (t) + 



k'k"k 



_ fc X! 5 k+k;k'+k" 5 ±+*^'+^"(-Sak)* (1) (*) 



ft 

cr'o-"5- 

k'k"k 



where, in order to lighten the notation, we dropped the two spin indexes a and a in the 
four-particle kernel function F defined in Eq. (1331) for they are already univocally deter- 
mined once chosen the others (i.e. o' and a") as soon as their selection rule (8 a+ a-y+a") is 
applied. Moreover, the H/2 has been reabsorbed in the Coulomb nonlinear coefficients V xx 
and F aa (t-f). A detail micro scop ic account for the mean-field V XX1 for the F's and their 



selection rules are considered in 



29 



301 ] . For the range of k-space of interest, i.e. |k| -C ^- 
(much lower than the inverse of the exciton Bohr radius) they are largely independent on the 
center of mass wave vectors. While and F ±=t (t — if) (i.e. co-circularly polarized waves) 

14 



conserve the polarizations, F ±T (t — t') and F T± (t — t') (counter-circular polarization) give 
rise to a mixing between the two circularly polarizations. The physical origin of the three 
terms in Eq. (1341) can be easily understood: the first is the Hartee-Fock or mean-field 
term representing the first order treatment in the Coulomb interaction between excitons, 
the second term is a pure biexciton (four-particle correlations) contribution. This coherent 
memory may be thought as a non-Markovian process involving the two-particle (excitons) 
states interacting with a bath of four-particle correlations |2j]. Equation (134"1) even if formally 
similar to that of Ref. [2J, represents its extension including polaritonic effects due to the 
presence of the cavity . It has been possible thanks to the inclusion of the dynamics of the 
cavity modes whereas in Ref. 2| the electromagnetic field entered as a parameter only. For- 



mer analogous extensions have been obtained within a semiclassical model [19j, [29|, [3Cj . The 
strong exciton-photon coupling does not modify the memory kernel because four-particle 
correlations do not couple directly to cavity photons. As pointed out clearly in Ref. l^ . 
cavity effects alter the phase dynamics of two-particle states during collisions, indeed, the 
phase of two-particle states in SMCs oscillates with a frequency which is modified respect to 
that of excitons in bare QWs, thus producing a modification of the integral in Eq. ( 13"4"1) . In 
this way the exciton-photon coupling V n u affects the exciton-exciton collisions that govern 
the polariton amplification process. Ref. 19] considers the first (mean-field) and the second 
(four-particle correlation) terms in the particular case of cocircularly polarized waves, calling 
them without indexes as and F(t) respectively. In Fig. 1 they show ^(uj), the Fourier 
transform of F(t) plus the mean-field term V xx , 



= V xx -i / dtF{t)e lut . (36) 

J — oo 

Its imaginary part is responsible for the frequency dependent excitation induced dephasing, 
it reflects the density of the states of two-exciton pair coherences. Towards the negative 
detuning region the dispersive part Re(jF) increases whereas the absorptive part Im(jF) 
goes to zero. The former comprises the mean-field contribution effectively reduced by the 
four-particle contribution. Indeed, the figure shows the case with a binding energy of 13.5 
meV, it gives V^xn sat — 11.39 meV which clearly is an upper bound for Re(jF) for negative 
detuning. The contribution carried by F(t) determines an effective reduction of the mean- 
field interaction (through its imaginary part which adds up to V xx ) and an excitation induced 
dephasing. It has been shown |l9j that both effects depends on the sum of the energies of 
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the scattered polariton pairs. The third term in Eq. (1341) can be thought as a reminder 
of the mismatch between the picture of a biexciton as a composite pair of exciton. In the 
following we will set ~ 2r x . 

The other nonlinear source term in Eq. (ETC]) depends directly on the exciton wave function 
and reads 

£<^>* (1) £ C^f k " (a k ^\B n , /k ,,)V . (37) 

nk n'k',n"k" 

It represents a phase-space filling (PSF) contribution, due to the Pauli blocking of electrons. 
It can be developed as follows, 

C th'nk' k " = V n'k'(^ak | 6 {n , k/) . ink) - [B nak ,Bl a , kl } | lnW) = 

= K'fc"W'< ^*naq*n' CT '(q+r, h (k'-k))(l^k | 4 '.q+^k'-kH^k'^Q+^k | lnW) + 

I q 

^^q$n' CT '(q-^(k'-k))(l^k | 4',-q+^(k'-k)+ % k'^,-q+%k | lnW) I = 

q J 
I q 

E $ -q^ CT 'q 3 ^q^n» CT "q 5 j, (38) 

q ) 

the explicit expressions of the q's are given in jsil ]. 

Thus, the nonlinear dynamics of Eq. (120]) driven by Hi can be written 



£/ £ \+(3) _ li^Ufh. I a , \+( 3 ) _ 1 A - (B 

, \ -Dno-k / ft h ^ ' h °k+k;k'+k" \ D h<r\d 



(It 



k'k"k 



(a CTk ') (1) (5n».k") (1) v;* [E + E $ n; $ ^ $ ij • (39) 

q q 

We are interested in studying polaritonic effects in SMCs where the optical response involves 
mainly excitons belonging to the IS band with wave vectors close to normal incidence, i.e. 
|k| <C ^- (much lower than the inverse of the exciton Bohr radius). In this case the exciton 
relative wave functions are independent on spins as well as on the center of mass wave 

vector. They are such that Y^=-oo l $ q| 2 = i- e - $ q = VaJuS^W' ° x is the exciton 
Bohr radius. From now on whenever no excitonic level is specified the IS label is understood. 
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It yields 



it < > 



"(3) 



.V^k , * v+(3) _ * A 
-1 ^ \ <Vk / ^ 2^ "k-l, I. I 

k'k"k 



(4k)* (1) (^) (1) (^k») (1) 2V;*0 PSF , (40) 

where the overlap PSF has been calculated in the case of zero center of mass wave vector, 
namely 

n 

In SMCs a measured parameter is the so-called vacuum Rabi splitting V na k |6[ of the IS 
excitonic resonance, for the range of k-space of interest essentially constant. Defining V = 
V a = V a VA<p*(0) 

(41) 



*,oPSF 



v*o 



V _ p SF = 8 nat v 1 V 



v/A/>*(0)~ 7 A ' 2n sat 

where we have set n sat = (7/16)-(A/7ra%), called saturation density. 

In terms of the two circular polarizations the dynamics induced by Hi finally reads 



d_ 

dt 



(B ±k ) 



+(3) 



Y 



v+(3) 



= +i-r ( a ±k ) 



E<4i)* <1, «±3 



(42) 



where 

E<s ±s >* (1) Cl = ^ E w~.<4 i > ,,,) fc 1 *> w <4*.>« • («) 

"■sat 

k k'k"k 

The same lines of argument can be followed for computing the Coulomb-induced interactions 

qCOUL(2) 

We are lead to introduce the saturation density for two main reasons. The most obvious 
is our interest to refer this work to the literature where n sat is extensively used 

BE" 
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32]. 



The other most interesting reason is that we can directly compute this quantity. Indeed, 
the equation of motion for the exciton operator reads 



< B ±k r< 3 > + ^ < a ±k )+< 3 > - 1 yj (B^r m Q^ 



V 



.\(«/6.,.\W 



-I 2 7S^° PSF E W*- +t »<B ±t )- l "<a ±k ,)'"<B ±k »> 

^ ' k'k"k 

Leaving apart the discrepancy between the order in the DCTS we can compute the so-called 
oscillator strength (OS), defined as what multiplies the photon expectation values (a± k =o}> 

2 



V 

OS =4— I I 



A(j)*(0) 



-O 



PSF 



(B ±0 y {1) (B ±0 ) 



(1) 



(44) 
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The saturation density may be defined as the exciton density that makes the oscillator 
strength to be zero. We obtain 

\VA<f)*(0) J ™l 16 

Eventually, the lowest order (x^) nonlinear optical response in SMCs are described by 
the following set of coupled equations: 

j t ( a ±k ) +,3) = -at ( a ±k ) +<3 » + ij ( « ±k ) +(3 » + t c ^f , (46) 
| ( s ±k >+*> = -a*, ( B±k }+«> + i\ < a ±k > +<3 » - 1 E^r'Mtt ■ (47) 

sic 

with E-c(5-t)* (1) i? (2) - = E-w(5-w)* (1) Q COU - (2) + Ew0+t)* (1) Q PSF -?, with the first of 
the two addenda originating from Coulomb interaction, Eq. fl35|) . whereas the second rep- 
resents the phase-space filling contribution written in Eqs. (l4"3{) . Starting from here, in 
the strong coupling case, it might be useful to transform the description into a polariton 
basis. The proper inclusion of dephasing/relaxation and the application of these equations 



to 



3arametric processes, in the strong coupling regime, is described in another paper of ours 



23]. 



Equations fj4"6j) and (|47|) is exact to the third order in the exciting field. While a sys- 
tematic treatment of higher-order optical nonlinearities would require an extension of the 
equations of motions (see e.g. Appendix), a restricted class of higher-order effects can be 
obtained from solving equations fj4"6]) and f|4T|) self-consistently up to arbitrary order as it 
is usually employed in standard nonlinear optics. This can be simply accomplished by re- 
placing, in the nonlinear sources, the linear excitonic polarization and light fields with the 



total fields 



19|, |22j. Multiple-scattering processes are expected to be very effective in 



cavity-embedded QWs due to multiple reflections at the Bragg mirrors. 

V. PARAMETRIC PHOTOLUMINESCENCE: TOWARDS SEMICONDUCTOR 
QUANTUM OPTICS 

Entanglement is one of the key features of quantum information and communication 



technology [4] . Parametric down-conversion is the most frequently used method to generate 
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highly entangled pairs of photons for quantum-optics applications, such as quantum cryp- 
tography and quantum teleportation. This optical nonlinear process consists of the 
scattering of two polaritons generated by a coherent pump beam into two final polariton 
modes. The total energy and momentum of the final pairs equal that of pump polariton 
pairs. The scattering can be spontaneous (parametric emission) or stimulated by a probe 
beam resonantly exciting one of the two final polariton modes. In 2005 an experiment 
probing quantum correlations of (parametrically emitted) cavity polaritons by exploiting 



quantum complementarity has been proposed and realized [ll|]. The most common set-up 
for parametric emission is the one where a single coherent pump feed resonantly excites the 
structure at a given energy and wave vector, k p . Within the DCTS we shall employ Eqs. 
( flTj) . (TT81) and Eq. (IT9"|) in operatorial form, provided all the equations to become fully 
significant as soon as the expectation value quantities we shall work out would lie within 



the consistent perturbative DCTS order we set from the beginning [26] . In order to be more 
specific we shall derive explicitly the case of input light beams activating only the IS" exciton 
sector with all the same circularly (e.g. a + ) polarization, thus excluding the coherent excita- 
tion of bound two-pair coherences (biexciton) mainly responsible for polarization-mixing Q . 
Equations involving polariton pairs with opposite polarization can be derived in complete 
analogy following the same steps. Starting from the Heisenberg equations for the exciton 
and photon operators and keeping only terms providing lowest order nonlinear response (in 
the input light field) we obtain, 

d V* * E 

—a k = -iu c k a k + i-^B k + t c — - , (48) 
at n n 

—B k = -iuj k B k + i— a k + 
at n 

k,k',a 

-~ ^ 2/3 - uj lk - u lk )(lk | B k | 2(3)X 1 ' k>0 X 0>2p . (49) 

fe/3 

In the following we will assume that the pump polaritons driven by a quite strong coherent 
input field consists of a classical (C-number) field. This approximation is in close resemblance 
to the two approximations performed in the previous section (linearization of fluctuations 
and coherent nonlinear processes). We shall show that under this approximation, we may 
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perform the same manipulations ending up to a set of coupled equations analogous to Eqs. 
fj46|) and (jUJ). In addition, having a precise set-up chosen, we will be able to specialize our 
equations and give an explicit account of the parametric contributions as well as the shifts 
the lowest order nonlinear dynamics provides. We shall retain only those terms containing 
the semiclassical pump amplitude at k p twice, thus focusing on the "direct" pump-induced 
nonlinear parametric scattering processes. It reads 



d - V 

— 5 ±k = -^ kJ B ±k + i— a± k - (50) 

~~fr~ \+k,k'+k''^l±k,0^0,l±k''>0;l±k'(^k'',k p 4',k p + <5 k ,kjA",k p + <%, k /k',kJ ~ 

Sa k,k',k" 

~J L ^ 2k /3 _ w ik ~~ w ik)(l°"k | £? ±k | 2apkp)X 1 ^oXo,2ff iL (8 kl 3,2kp + h,k p $k g ,k+k P ) > 

where we have already manipulated the phase- space filling matrix element. Here in brack- 
ets the first addendum of each line would be responsible for the parametric contribution, 
whereas the others will give the shifts. It is understood, from now on, that the pump-driven 
terms (e.g. the X and Y at k p ) are C-numbers coherent amplitudes like the semiclassi- 
cal electromagnetic pump field, we will make such distinction in marking with a "hat" the 
operators only. 

We need some care in manipulating the Coulomb- induced terms, the last line. Written 
explicitly it is 



Coul 



~\ ~ w ik ~ ^ik)(l^k | B ±k I 2apkp)X ia ^ X ,2a p u p + 

ali.,apkp 

~\ ^ 2k /3 ~~ w ik ~~ ^ik)(lc r k | 5 ±k | 2o- /3 k /3 )X 1(7pkpi oXo i 2 CTk+kj) (k+k P ) (51) 

As for the term containing Xo t 2k p , we are facing a C-number which gives no problem in 
performing the very same procedure of the previous chapter. As for the other we would 
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exploit the formal biexciton solution 

X 0;2(fc+M (t) = j dt'e-^^^ 1 - (vkpMh + k p ) | B[ p I lk)X 0tlk Y 0Akp + 

V k (2(k + k p ) | B\ I lk p )X 0)lk Y 0>lk ^ , (52) 

where, for the sake of consistence, we are neglecting X .2( k+kp )(— oo) because the biexciton, 
within the present approximations, is always generated by an operator at k times a classical 
amplitude at k p which is always zero before the electromagnetic impulse arrived. Moreover, 
an analogous identity such that of Eq. (l2~5|) is valid in the present context, namely 



d 
dt 



Xo^Xo^e-^+^^-^j = (53) 

With these tools at hand we are able to perform step by step the manipulations of the 
previous section for all the quantities at play. The final result reads 

d * h y A 

i V ( * 

~Ti~n \ 1±k *>° °> 1±k 3> O' 1 * 1 '*' ^l±k p ,0^0,l±kp^0,l±k + ^l±k p ,0^0,l±k^0,l±kp 

-^ 1±ki ,o(t)|K x Xo, 1±kp (t)X , 1±kp (t) - i j dt'F^it - O^oWOAWO 

-2^X lCTkpkpi o(t)|K x Xo, 1±k (t)X , 1±kp (t) - z y" rft / F ±± (t - 0*0,1^(0*0,^(0 J (54) 

where kj = 2k p — k, and again and F ±=t (t — t') have reabsorbed the 1/2 originating 
from Eq. (153]) . In the specific case under analysis we are considering co-circularly polarized 
waves and the mean field term, as well as the the kernel function F(t) can be found in 



Refs. 



29 



30J. 



Eventually, the lowest order (x^) nonlinear optical response in SMCs is given by the 
following set of coupled equations where, in the same spirit of the final remark in the 
previous section, we account for multiple scattering simply by replacing the linear excitonic 
polarization and light fields with the total fields: 

— a±k = -zu; k a ±k + %— B± k + t c —— 
dt n n 

— B± k = -iuj^B^ + s ±k + i— a ±k - ^-R±k , (55) 
21 



where = (R%£ + R^ k ) 



sat v - 1 

R±k = - — -B±k p a ±kp 5 ±k 

"-sat 



R±* = Bl ki (t) (v^B ±kp (t)B ±kp (t) - 

-if dt'F^{t-t')B ± ^{t')B ±K {t')\. (56) 
The pump induced renormalization of the exciton dispersion gives a frequency shift 

./ v 

3±k = -i\ (5ik p a ±k p -B±k + B^ p B ±kp a ±]l ) + 

V— , 

^-^--Dikp^ikp^ik - 

-2^* kp (t) jT dt'F ±± (f - 0^(04^(0) • (57) 

Equations, fl55l) are the main result of this paper. They can be considered the starting 
point for the microscopic description of quantum optical effects in SMCs. These equa- 
tions extend the usual semiclassical description of Coulomb interaction effects, in terms of 
a mean-field term plus a genuine non-instantaneous four-particle correlation, to quantum 
optical effects. Analogous equations can be obtained starting from an effective Hamiltonian 

n 

describing excitons as interacting bosons [10]. The resulting equations (usually developed 
in a polariton basis) do not include correlation effects beyond Hartree-Fock. Moreover the 
interaction terms due to phase space filling differs from those obtaind within the present 
approach not based on an effective Hamiltonian. 

Only the many-body electronic Hamiltonian, the intracavity-photon Hamiltonian and 
the Hamiltonian describing their mutual interaction have been taken into account. Losses 



through mirrors, decoherence and noise due to environment interactions as we 
tions of this theoretical framework will be addressed in another paper of ours 



1 as applica- 



23|]. 



VI. CONCLUSION 



In this paper we set a dynamics controlled truncation scheme approach to nonlinear 
optical processes in cavity embedded semiconductor QWs without any assumption on the 
quantum statistics of the excitons involved. This approach represents the starting point 
for the microscopic analysis to quantum optics experiments in the strong coupling regime. 
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We presented a systematic theory of Coulomb-induced correlation effects in the nonlinear 
optical processes in SMCs. We end up with dynamical equations for exciton and photon 
operators which extend the usual semiclassical description of Coulomb interaction effects, 
in terms of a mean-field term plus a genuine non-instantaneous four-particle correlation, to 
quantum optical effects. The proper inclusion of the detrimental environment interactions 
as well as applications of the present theoretical scheme will be presented in another paper 



of ours 



23|]. 
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APPENDIX A: THE EQUATION OF MOTION AT ANY ORDER 



The equation of motion for the operators in (flOl) . under the Hamiltonian H = H e + H c + 
Hi + H p reads: 



d 



3=1 



8=1 



Na;MP 



kv 



Ek- 



6(m - 2)(m/i | a[ \ (m - —Y nX -( m -i)u + 



Ey 



6(n - 2)((n - | % | nA)— Y" (n _i )l/;mM 
-9(m - l)(m/i | as | (m + l)^)-^K n Ai(m+i)i/ - 

I \ \ ^ \ 

"fc I (n+l)i/;m/x 



e(n- l)((n + l)z/ | al | nA)-fF ( 



+ 



+ij mil e(M - 1) ^ (M/3 | ^ | (M + 



Afa;(Af+l)<5^nA,0 — 



7(5 



U n>1 e(N V ^ N + ^ I B lx I Na)X (N+1)v . tM/3 Y . m/ , 



Un, 5 x , Q(N -2)J2 VMN - 1)77 I 4s I Na)X { 



N~l)r);Mf3Y lk . mlM + 



nkr) 



+^<W><W>e(M -2)J2 Vnk(MP | B\- k | (M - l)5)X^ ;(A/ _ 1)5 f„A ; iS + 



nk5 



nku 



-T$N,o8 a ,o®(n - 2)^Fa((n - 1)7 I a k \ nX)X ln - k . M/3 Y {n _ lh . mfl 

nk'y 

-^N,i®{n - l)V a ^((n+ 1)7 I 4 a I r2A)X . M/3 Y" (n+1)7;mM + 
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nk vo 



V*- k (6(M - l)9(m - 2)(M/3 \ B nk \ (M + 1)6) 

(m/i | a\ | (m - l)u)X Na .^M+i)sY n x;(m-i)u - 
-Q(N - 2)6(n - l)((iV - 1)5 | B nl \ Na) 
((n + l)z/ | at | n\)X {N _ l)5 . M pY {n+l)u . mi ^ + 

^(e(iV - l)9(n - 2)({N + 1)5 I B\- k \ Na) 

((n — l)z/ | | nA)X( A r +1 ) (5;M/ gy( n _ 1 )j y;mM — 
-6(M - 2)0(m - 1)(M[3 | | (M - 1)5) 



(m/i | % | (m+ l)z/)X A r a;(M „ 1 ) (5 '5 > nA;(m+1)l 



(Al) 



Here Q(x) is the Heaviside function equal to 1 for positive argument and zero otherwise. 
APPENDIX B: N EH PAIR AGGREGATES 

We start from the usual model for the electronic Hamiltonian of a direct two-band semi- 
conductor 15j,|24|. It is obtained from the many-body Hamiltonian of the interacting electron 



system in a lattice, keeping explicitly only those terms in the Coulomb interaction preserving 
the number of electrons in a given band and can be expressed as 



H e = Hq + Vcoul 



(Bl) 



It comprises the single-particle Hamiltonian terms for electrons in conduction band and 
holes in valence band (here k = (k, a) and c CT) k (flV.k) annihilates an electron (a hole)) : 

H = J2 E c,kc[ck + E h,kd\dk , (B2) 

k k 

and the Coulomb interaction term of three contributions, the two repulsive electron-electron 
(e-e) and hole-hole (h-h) terms and the attractive (e-h) one: 



Vc 



Coul 



lj2 V 1 Yl Ck+q^k'-q^k'^k + ^ 4,k+q4',k'-q^',k'4k ~ 



q^O <x,k,<r',k' q^O <7,k,<r',k' 

-/] V 9 ^2 4,k+q4',k'-q^',k'C CT ,k • 



(B3) 



q^O a-,k,o-',k' 

A many-body interacting state is usually very different from a product state, however a 
common way to express the former is by a superposition of uncorrelated product states. 
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The physical picture that arises out of it expresses the dressing the interaction performs 
over a set of noninteracting particles. The general many-body Schrodinger equation for this 
Coulomb-correlated system is 

H e | = (H + Vcoui) | *> = E | *) , (B4) 

with | the global interacting many-body state of the whole Fock space and E its corre- 
sponding energy. The system Hamiltonian commutes with the total-number operators for 
electron and holes, i.e. N e = "^2 k c] e Ck and Nh = Sfc^l^fc- Therefore the state | may be 
build up corresponding on a given number of electrons and of holes. Moreover, because we 
shall consider the case of intrinsic semiconductors materials where N e = = N, the good 
quantum number for the Schrodinger equation (1B4I) is the total number of electron-hole 
pairs N, explicitly 

H e | Na) = E Na | Na) , (B5) 

where a is the whole set of proper quantum numbers needed to specify univocally the many- 
body state. 

For any given number N of electron-hole pairs, the product-state set, built up from the 
single-particle states {| Na)} eigenstates of the noninteracting carrier Hamiltonian H , is a 
natural complete basis of the N-pair sector of the global Fock space: 

H | Na) = e Na | Na) , (B6) 

where N identifies the N-pair subspace and a is a compact form for all the single particle 
indexes, i.e. a = j el , j e2 , j eN ; j hl , j h2 , j h N- Indeed 

N 

I N, a) = ®^ =l c)J) hn | 0, 0) and e Na = ^(e jen + e jhn ) . (B7) 

n=l 

Being a complete orthonormal basis for the N-pair subspace we may expand the many-body 
state | N, a) over it, it yields 

\Na) =J2 U a a I Na ) ■ ( B8 ) 

a 

It is only a matter of calculation to show that U^ a is nothing but the envelope function of 
the N-pair aggregate, solution of the corresponding secular equation. Indeed the eigenvalue 
problem ( 1B5I) is transformed into: 

J2(( Na I ^ | Na') - E Na 5 a , a ,)U$ a = . (B9) 

a' 
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Namely N=l leads to the exciton secular equation, whereas N=2 represents the biexciton 
(two pairs) Coulomb problem. 

In order to be clearer we shall propose in details the N=l exciton calculation. We 
shall work in the direct lattice r <-> r, (the former is a continuous variable whereas the 

n 

latter is a point in the 3D lattice). Using the general mapping [33j Yl r - (V^o) / d 3r i 
5(r — r') = (8 ri ,T S /vo), and a\. = (^(r^/y/vo), here vq is the unit cell volume and for 
simplicity the spin selection rules for the optically active states has been already taken into 
account, (1B9|) reads 

J2 ( r h | H e | r' e , v' h ) - JWr.^ ) f/ a Q (r' e , A) = , (BIO) 



r' r' 



with 



f) 2 h 2 e 2 



(r„r, | H. | r>A) = | - — Vj. - —V th - -j—j + V M | . (BU) 

here V(r e , r^) represents all the additional potential, e.g. those of the heterostructures 
or those of disorder effects, V(r e , r^) = V e (z e ) + V h {zh). Typically, the energy difference 
between the lowest QW subband level and the first excited one (at least a few meV) is much 
larger than the Coulomb interaction between particles (a few meV). As a consequence, 
at least at low temperatures, particles are confined at the lowest quantization level and 
the (possible) distorsion of the wave function due to the Coulomb-activated admixture of 
different subbands can be safely neglected. In some extent, then, the particle wave function 
dependence along the growth (say z) direction can be factorized out and the dynamics 
becomes essentially two-dimensional. However, a purely 2D approximation for excitons 
would miss important effects of the geometrical QWs parameters on the binding energy 
and would not be able to account for the interaction with a 3D continuum environment 
of surrounding modes (e.g. acoustic phonon modes in heterostructures with alloy lattice 
constant in close proximity [34]). In addition in QWs, light and heavy holes in valence band 
are split off in energy. Assuming that this splitting is much larger than kinetic energies of 
all the involved particles and, as well, much larger than the interaction between them, we 
shall consider only heavy hole states as occupied. 

In Eq. fIBlOp the 3D Coulomb interaction prevents form factorizing into (free) in-plane 
and confined directions. Nevertheless if we assume that the quantization energy along z is 
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much larger than the Coulomb energy, at leading order we can factorize out the ^-dependence 

h 2 d 2 , h 2 d 2 



+ V e (z e )-—— T + V h (z h ) \U a (r e ,r h ) 



2m e dz 2 6 2rrih dz 2 

E z U a (r e ,r h ). (B12) 

It means we are solving our secular equation with solutions built up as linear combination 
of F£ n a (r|, rj l )c nc (z e )v nv (zh), with r = (r\z). Equation flB12j) expresses the lack of 
translational symmetry along the growth z direction, thus single particle states experience 
confinement and two additional QW subband quantum numbers n v , n c (for valence and 
conduction states respectively) appear. We still leave a as a reminder that new possible 
indexes could still arise in due course. 



Projecting Eq. ( 1B10I) on these confined states we end up with an effective Schrodinger 
equation in the plane 



= (E a - E z nc - E z nv )F: M (rl rj) , (B13) 

with 

U nc ,n v ;n' c ,n' v (\ r l ~ r hl) = / dz e / dz h & -.C nc (z e )c n > c {z e )v Uv (z h )v K (z h ) . 

J J e r ^\vl-v{\ 2 + {z e -z h ) 2 

(B14) 

For what already stated, we shall consider only the lowest confined subband levels, then the 
resulting effective in-plane secular equation becomes 

2m e r e 2mh r h J J , . /,J| _ r JI p + ^ _ ^ J 

{E a -E z c -E z v )FZ{viv{), 



with the product exciton envelope function U a (r e ,r h ) = F"(ri, r\)c(z e )v(zh)- Equation 
(IB15P is solvable by separation of variables once we employ a coordinate transformation into 
center of mass (CM) R = (m e rl + mhv\) / {m e + rrih) and relative p = (r| — rj[) exciton 
coordinates. It reads 



Vk- TT-^l ~ U(p) \F:(R,p) = EF:(R,p) , (B16) 
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with a solution we can arrange as F°(R, p) = £t ^ W^(p) the latter solution of the relative 
hydrogen- like 2D problem. 

Eventually, in real-space representation, we have our exciton wave function with total 
in-plane CM wave vector K (A is the in-plane quantization surface in the free directions) 
which reads 

I n(xK > = 7^ E elk ' R WUp)c(zMz h )al e air h I 0) , (B17) 

r e ,r h 

being a^ VT (a c / V;T ) creation (annihilation) operator of the conduction- or valence- band elec- 
tron in the Wannier representation and r e /h = (rj!^, z e /h) are to be considered coordinates 
of the direct lattice, | 0) is the crystal ground state. 

When e.g. exploring the exciton-phonon interaction, it is useful to express exciton states 
in reciprocal space. With the usual transformation to Bloch representation, (N = v A L is 
the number of unit cells and L is the quantization dimension along the confined direction, 
v = c,v), 

au(r, z) = ^ E e ' k ' rn ^(k,M , ( B1 8) 

one obtains: 

I naK ^ = E 6 *^' j d P J dz -\ ^^(p)c(^)^(^)e-^- ( ^ k+J ' ek ' ) e-^^ e e - J ^^ 

k z ,k z 

a l,(Kk z ) a ^',K)\°) ■ ( B19 ) 

In order to end up with a form as much as possible in analogy with its bulk counterpart we 
shall define CM and relative coordinates even in the reciprocal lattice: 

K = k - k' f k = k r + rj e K 



(B20) 

k r = r] h k + i] e k' I k' = k r - i] h K 



It becomes 



K) = <W-k' E i^jj dpW na (p)e- lp ^) ( -j= J dz e c(z e )e- ik > 

-±= J dz h v(z h )e- ik '^ J< (kr+ ^K 1 fc,) a «.(kr-%K 1 fc;)|0> . (B21) 



Thus 



naK) = J2J2 $ n,,k,<^<( kr+ , e K, fe ,)^(k r -^K, fci )| > > (B22) 
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or in the electron-hole picture (a Vj k = (f_ k and — k \ el = k \ hole ) 



no 



K> = ££* 



K n e ii h 

na^ r a k z u k' z L '(\ lr +ri e K,k z ) u '(-kr+rih'K-k' z ) I 



) d \-k r +il h -K,-k')\°) ' 



(B23) 



with the relations in Eq. flB20j) changed accordingly. 
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